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I. INTRODUCTION 

In this paper a quantitative understanding of Kaon-few-nucleon quasi-bound states is attempted. In recent years, 
the existence of such states has been vividly discussed. It was initiated by the KEK finding of peaks in the nucleon 
spectra of K~ absorption in 4 He, [HQ. Additional evidence was given by the FINUDA measurement of the invariant 
mass distribution of the Ap produced in K~ absorption by light nuclei [3|. The existence of such bound states have 
been expected as the kaon- nucleon and the kaon- nucleus interactions have been known to be strongly attractive, (4j. 
This is now firmly confirmed on the basis of kaonic atom data [5J. However, the KEK and FINUDA experiments 
indicate unexpectedly strong bindings of the order of 100, 150 MeV in the lightest KNN, KNNN systems. These 
experiments require further confirmation. Also, the interpretation of the observed peaks has been disputed in Refs. 
@, [3] while the initial interpretation is defended in Ref. 

Calculations indicate that such states are expected, albeit these might be very broad and difficult to detect. The first 
calculations performed by Akaishi and Yamazaki in Ref. Q were followed by several subsequent publications. These 
calculations exploited essentially the S - wave resonant attraction related to the A(1405) state. With an optical model 
type of approach it was shown that the K-meson optical potential at the center of small nuclei may be as strong as 
500 MeV generating very strong binding of the meson and a strong contraction of the few-nucleon systems. However, 
to reproduce the KEK data, these calculations involved some relaxation of the NN repulsion at short distances which 
would allow the existence of strongly bound and very dense systems. These calculations raise the important question 
on how to implement a realistic short range NN repulsion in the kaonic systems. 

Another open question is related to the strength and range of KN interactions. Any mathematical description of 
few body systems requires knowledge of NN and KN off-shell scattering amplitudes. Those related to NN interactions 
are controlled fairly well in terms of modern NN potentials. For a bound K-meson the amplitudes needed involve the 
subthreshold energy region 

JkN = Jkn( — Eb — -Erccoil), (1) 

where Eb is the KN separation energy and -E r0C oii the recoil energy of the KN pair relative to the rest of the system. 
If the separation energy is as large as 100 MeV, meson momenta become 250 MeV/c and E xeco i\ may be as large as 
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40 MeV. The energies of interest for (—Eb — -^recoil) ar e then located well below the A(1405) state. The amplitudes 
there are strongly attractive and so when used in a standard optical potential approach may support very strong 
bindings. One problem that arises at this stage is of a technical character. As these amplitudes are energy dependent, 
it is hard to account for that in the optical model approach. There exists another, more serious, problem which is 
common to all approaches. As the energies involved are far away from the physical region tested in KN scattering, the 
uncertainties in the KN scattering amplitudes are sizable. For instance, if the A(1405) is a KN bound state, then the 
amplitude far below the resonance is given not only by the position of the singularity but to a greater extent by the 
Born term, which indicates a strong dependence on the uncertain interaction range r D . An old multi-channel potential 
model of Ref. [ll| indicates that the available scattering data do not allow one to fix the precise value of r Q . This 
unfortunate situation is still actual. The r a is expected to be close to the inverse vector meson mass. However, even 
though a change of 20% in r a would not affect the scattering region, it results in a 30% change of Jkn m the deep 
subthreshold region. The corresponding uncertainty in the binding energy then amounts to ~ 30 MeV. As indicated 
by few body calculations of Ref. [lfj this problem strongly affects the outcome. The uncertainties of Jkn require 
further coherent experimental and theoretical studies of KN and K-few-N interactions. It becomes one of the most 
important purposes of the K meson physics. 

On the other hand, there is one consequence of Eq.([T]) which is model independent. If the binding and recoil 
are so large the fxN amplitudes involve the energies below the thresholds of meson- hyperon decay channels. As a 
consequence the dominant decay modes are blocked and the lifetimes of nuclear K meson systems are determined only 
by multi-nucleon captures. This leads to the expectation that such states may live long enough to be detectable, (4J. 

There exists several calculations of KNN binding energies. These states are named K^pp although in reality they 
correspond to isospin /jvjv = 1 and total isospin I tot = 1/2. The first prediction by Akaishi and Yamazaki lead to 
(E B ,T) = (48,60) MeV 0. With a similar, molecular type method, Dote and Weise [U[ obtain E B , < 50 MeV and 
indicate a strong dependence of this result on the short range NN repulsion. On the other hand, the recent three 
body calculations based on Faddeev or AGS methods yield larger bindings. Thus Schevchenko et al. [IH obtain 
(E B , T) = (55 - 70, 95 - 110) MeV while Ikeda and Sato [3 calculate (Eb, L) = (~ 80, - 73) MeV. Later in the text 
we show that the discrepancy of these two groups of results is due to: different A(1405) properties, explicit description 
of the multiple scattering in decay channels and possibly to an incompatible treatment of the NN repulsion. 

There are two new elements introduced in this paper. First, the P-wave interactions due to £(1385) have been 
indicated as a possible source of the strong binding. Here, these are introduced explicitly. Second, the stress is put 
on the strong KN spacial correlations induced by the S and P wave resonances. 

Leaving aside the interpretation of the peaks attributed to bound KNN and KNNN systems the essential theoretical 
questions are: 

(a) What is the binding mechanism? 

(b) Are the technical questions under control? 

(c) Can the widths be narrow? 

This paper attempts an answer to these questions and the following results are obtained : 

1) To account properly for the KN force range, short range KN correlations and the NN repulsion, a two step 
calculation is performed. First a wave function involving strongly correlated K-N subsystems is found in a fixed 
nucleon approximation. This step also allows one to find potentials due to the K meson which tend to contract 
the nucleons. Next, these correlated wave functions and contracting potentials are used as the input in variational 
calculations for the K-few nucleon binding. In the KNN case the binding energy and width are found by solving the 
Schrodinger equation. 

2) While the dominant mechanism of attraction is related to the A(1405) state, it is found that another resonant 
state, the £(1385), contributes significantly to the structure of the bound states but much less to the binding in KNN 
and K-few-N systems. In addition the £(1385) generates new branches of nuclear states that could not be generated 
by the A(1405) alone. 

3) The binding energy is determined to a large extent by the attraction and the repulsive core in NN interactions. 
With the Argonne NN potential [H[ one obtains the lowest state of KNN bound by about 40-80 MeV and a KNNN state 
bound by about 90-150 MeV. Moderate dependence on the KN interactions is found, provided these are constrained 
by the shape of A(1405) and the value of the KN scattering length. However, the position of A(1405) itself is not well 
known and this becomes the source of a large uncertainty. The effect of £(1385) on the binding energy is limited, 
In the states bound via A(1405) it adds some 5-10 MeV contribution to the KNN binding and 10-20 MeV to KNNN 
binding. In this sense the suggestions of Ref. [l(| are not fully supported. However, the effect of £(1385) on the 
space structure of deeply bound kaonic states is strong. The £(1385) is formed in peripheral regions and it absorbs a 
large fraction of the total K meson binding. In consequence the radii of these systems are fairly large and the nucleon 
densities are comparable to those met in the 4 He nuclei. 

4) The problem of uncertainties related to the large recoil momenta entering Eq.([T]) is only partly removed. Large 
kaon momenta are hidden inside the resonant structures. In principle these may be kept under control with the help 
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of other experiments. In practice it is not the case. The other sector of large momenta, due to the strong binding, is 
partly screened by the short range NN repulsion. The main consequence is a strong dependence of the meson binding 
energies on the position of the A(1405) resonance. In principle the shape of A(1405) is tested by the invariant mass 
distribution in the decay E-7T channel. In practice it is not so as the relevant energy region is located close to the T,n 
threshold. In this region the theoretical and experimental uncertainties are large. 

5) These states are very broad if the binding energies are less than 100 MeV. For stronger bindings, which are 
possible under the current values of the KN parameters the main mesonic decay modes may be closed. The widths 
for non-mesonic modes are hard to calculate and extrapolations from the emulsion data are not very reliable. New 
experiments are needed. 

A simple physical picture emerges from this approach. The mesons are strongly correlated to slowly moving 
nucleons. The correlations are of the A(1405) type at large densities, and of the £(1385) type in the peripheries. Each 
K,N pair has a good chance to stay also in the E,7r form. The structure is rather loose as sizable fractions of the 
binding energies are hidden in the short ranged correlations. 



II. THE KNN BOUND STATE 



This section presents an introduction to the method used in this work. Several steps describe the increasing degree 
of precision and also the increasing level of technical complications: 

• At first the KNN levels are found within the fixed nucleon approximation with a simple S wave KN interaction. 

• The nucleon degrees of freedom and NN interactions are introduced and a related Schrodinger equation is solved. 

• The method is extended to multiple channel situations. 

• Both S and P wave KN interactions are allowed. 

Consider scattering of a light meson on two identical, heavy nucleons. To begin with, the nucleons are fixed at 
coordinates Xj(« = 1,2) and the wave function is assumed to be in the form 



*(x, Xi,x 2 )=Xi^( x ; Xi,x 2 ) XNAf(xi,x 2 ), (2) 

where x is the meson coordinate. The notation is simplified and some possible indices are suppressed. The meson 
wave function \K is given by the solution of the multiple scattering equation 

/exp lip I x _ y 1 1 
dy — — ; j— U K N(y,Xi) Xfir(y,xi,x 2 ) (3) 
4tt I x-y | 

obtained with fixed positions of the nucleons. An equation of similar structure with a zero range meson-nucleon 
pseudo-potential U was used by Brueckner [13] to calculate the scattering length of a meson on two nucleons. For a 
high energy scattering it was extensively discussed by Foldy and Walecka, who used finite range separable interactions 
U [HI]. With such interactions equation ([3]) allows for semi-analytic solutions in the NN, and also in few nucleon 
cases. Here, the method is extended to the bound state problem. One looks for solutions of Eq. |S]) with no incident 
wave xk ( x )°- The momentum p becomes a complex eigenvalue p(xi) which determines the energy and width of the 
system for given nucleon positions Xj. 

Equation (|3|) is written in terms of the Klein-Gordon or Schrodinger propagator. The difference arises when the 
relation of energy and momentum is established. Reasons of simplicity, which will become clear later, favor the non- 
relativistic relation in the KN center of mass system. Thus, the interaction is presented as Ukn = Z^knVkn where 
fiKN is the reduced mass. Corrections for relativity may be introduced at a later stage. The potential Vkn for an S 
wave interaction is chosen in a separable form 

Vfav(x - x ; ,x' - xi) = A u(x - xi) v(x' - xi), (4) 

where v is a form- factor and A is a strength parameter. The eigenvalue equation is now reduced to 

/ x , , f, exp[ip(x 1 ,x 2 ) I x-y |] f , , , 

Xk(x,Xi,x 2 ) + £* A J dy — j - _ ■ v(y - x ; ) J dy' v{y - x s ) x K (y , xi, x 2 ) = 0. (5) 

Equation (f5J) becomes a matrix equation for wave amplitudes ipi defined at each scattcrcr i by 

tpi=X J dx v(x- xi) xk(x,x 1 ,x 2 ). (6) 
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To find the equations for ipi one introduces the off-shell KN scattering matrices / and matrix elements of the propagator 

n < \ [a a i , cxp(ik | x — y |) 

G«: J (x i ,Xj) = y dydxw(x-Xi) ^ - _ - u(y-Xj). (7) 

The diagonal value, Gj^ = G, determines the meson nucleon scattering matrix t by the well known (see e.g. Ref. (l8| ) 
relation 

= (1 + A G)- 1 A (8) 

and this yields the full off-shell scattering amplitude / 

f(k,E,k')=v(k)t(E)v(k'). (9) 

Here, k, k' are the initial and final momenta while the form-factor v{k) is given by the Fourier transform of v(r). The 
Yamaguchi form v(k) = 1/(1 + k 2 / k 2 ) with a free parameter k will be used in this paper. At zero momenta and at the 
threshold this choice normalizes / (and t) to the scattering length. Unfortunately for historical reasons the standard 
convention in the K-N system is to define the scattering length by 

a + ib = -f(k = 0, E = 0, k' = 0) = F(0, 0, 0) (10) 

and the capital F will be used in several places to comply with the standard KN parameters. 

In order to cast Eq. © into a standard multiple scattering equation for ipi one carries out the following three steps: 
1) Integrates Eq. (JSJ over the i-th form-factor u(x — Xi). 2) Selects the i-th term from the R.H. side. 3) Multiplies 
Eq. ||SJ) by (1 + AG) -1 . In this way the kernel of the multiple scattering equation can be expressed in terms of 
scattering amplitudes t% at each nucleon i and propagators describing the passage from the nucleon i to the other 
nucleon j . One now arrives at a set of linear equations 

ipi + Sj^j tj G it j tpj = 0, (11) 

which may be solved numerically. For the Yamaguchi form-factors, propagators Gij allow analytic expressions 

1 k 2 + k 2 

Gi 2{r, k) — — v(k) 2 [exp(ifcr) — exp(— nr) — r — exp(— nr)] = G(r,k), (12) 

r 2k 

where r = X2 — xi. For the sake of illustration, the KNN case is presented in some detail. The condition for a bound 
state with two amplitudes leads to a pair of equations 

V>i + 1 G V> 2 = 0, ip 2 + t G i/n = 0. (13) 

When the determinant 

D = 1 - (t G) 2 (14) 

is put to zero, the binding "momenta" p(r) may be obtained numerically. Two different solutions corresponding to 
1 + tG = or 1 — tG = may exist. The first solution is symmetric ip2 — ipi and describes the meson in the S wave 
state with respect to the NN center of mass. The second solution is antisymmetric ip2 = — ipi and describes a P wave 
solution. With the rank one separable interaction this latter solution does not exist in the full range of r. However, 
it arises with the more complicated rank two interactions discussed later. 

Eigenvalues corresponding to unstable quasi-bound states are obtained in the second quadrant of complex p(r) = 
p = Pr + ipi plane. In this quadrant the kernel 

2 2 

tG = f(p) [exp(-p 7 r) exp(ip R r) - exp(-Kr) (1 + r — )]/r (15) 

Ik 

is exponentially damped at large distances as required by the asymptotic form of the bound state wave function \k- 
At short distances G is regularized by the KN form-factor. 

If the scattering amplitude is dominated by a quasi-bound state, such as A(1405), the related pole dominates and 
in some energy region / ~ 7 2 j(E — E*), where 7 is a coupling constant and E* — E r — iT r /2 is a complex A(1405) 
binding energy. The full KNN binding energy, Vk, is given by the equation 1 + tG = which becomes 

V K {r)~E* - 7 2 G(r,p). (16) 
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Since Re G(r,p) close to the resonance is positive this solution offers binding stronger than the « 28 MeV binding 
in the A(1405). Asymptotically, for r — > oo one obtains Vk — * that is a kaon bound to a nucleon to form a 
A(1405). This type of asymptotic behavior occurs in all of the K-few-nucleon systems of practical interest. Hence, 
the separation energy is understood here as the separation of the K-N-N system into the N-A(1405) system. 

The limits r — ► in Eqs. (| 1511 6[) are regular. However, a joint limit of zero range KN interactions, K — * oo, 
and r — > is singular and the KNN system collapses. Therefore, some care is necessary when this limit is taken. 
Here, we stay within a phenomenological approach and the standard expectation that the range of KN interactions is 
determined by vector meson exchange. In equation (|12[) for G the range of interactions enters twice, first as a cutoff 
at small distances and second in terms of the form- factor v(k) 2 . We find in a numerical way that these two effects 
cancel and fG is very stable within the range 3 < K < 6 fm. As the KN interaction range is very short, but finite, 
the uncertainties related to the actual value of n arc additionally eliminated by the short range repulsion in the NN 
systems. This yields an important stability in the few-body calculations described here. 

With nucleons fixed at a distance r the eigen-value condition determines p(r) which in turn generates the potential 
Vk(t) contracting the NN system to a smaller radius. The form and strength of this potential depends on the form 
of the kinetic energy. In the KN CM. system the Ekn energy is given by \J M 2 + p 2 + \J m 2 + p 2 , where m, M are 
masses of the meson and nucleon. The same form is kept in the large nucleon mass M limit. In the few nucleon 
systems the non-relativistic form of the nucleon energy is used and the problem of large nucleon mass disappears. 

The meson propagator in Eq. ^ is chosen to make the multiple scattering equation (jlip equivalent to a differential 
equation 

[-A + 2/iKiv V KNi ] X = P(r? (17) 

i 

and the contracting potential becomes Vk = p(r) 2 /2^lkn ■ The advantage of this choice is discussed in the next 
section. 



A. Schrodinger equation 

The solution of the full KNN bound state problem is given by equation 

("^ ~ ^7 " ^7 + Vkn1 + VkN2 + Vnn ^ = E *- ( 18 ) 
Zm ZM ZM 

The wave function is assumed in the form — xk(x, Xi)XNN{r) as given in Eq. Multiplying Eq. fTS]) on the left 
by Xk and integrating over the meson coordinate x one obtains the Schrodinger equation for the NN wave function 



Xnn(v) = J dx XK(x,Xi)^(x,Xi) (19) 

in the form 

[E - V K (r) + Ax/2M + A 2 /2M - V NN ] X nn + AE kinXN N = 0. (20) 

where the last term AEkin is a correction to the kinetic energies. This correction is small due to the choice of the 
meson kinetic energies. In the Schrodinger equation (fl8|) it is given by the meson mass m. On the other hand, to 
determine the A(1405) properties and to solve the scattering equation ([3|) the reduced mass hkn is used. Due to this, 
the correction term AEkin is of the order of l/M. In addition, the meson wave function satisfies the relation 

A xX if =^A iX k, (21) 

i 

which may be obtained by partial integration over coordinate y in Eq. ©. In this way 

AE kin XNN = J dx XK ~3iXK ~^iXNN, (22) 

which is very small due to angular averaging and sign changes in the derivatives. In more detail this correction reduces 
to 

AE MK . -A J d( IT p^|z^g>_ ^ NN(r) , (23) 
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and is suppressed by the angular average over £ and at large r by the small overlap of G(£ — r) and G(£). The AE^in 
makes a contribution w 0.2 MeV to the binding energy. Such twice damped, small terms of similar type, arise also in 
more involved versions of this calculation. The AEkin is of the same order but is given by very lengthy formulas. Since 
it is very small in comparison to the dominant uncertainties in Vk it is dropped, leading to a significant simplification 
of the variational approach. 

As the next step, equation ([20]) is solved with an S-wave interaction based on the more realistic NN potential 
of Argonne [HI]. This solution is also compared to another, variational solution with the intention of checking the 
variational method used in heavier systems. The actual interaction used, in the notation of Ref. [ll|, has the form 

v(NN) = v EM (NN) + v v (NN) + v R (NN), (24) 

where the electromagnetic part v EM only includes the dominant term proportional to Fc(r) in Eq. (4) of Ref. [ill ], 
the OPE term is given by Eq. (18) of Ref. [l5[ and the phenomenological short range term v R from Eq. (20) 
with the parameters in Table II - again all from Ref. [lj|. This gives directly the S-wave T=l, S=0 interaction 
v(S — wave, T = 1, S = 0). However, in the T=0, S=l deuteron channel, the effect of the tensor interaction 
vt(T — 0, S — 1) on the central component v c (T = 0, S = 1) is incorporated by the closure approximation to give 

V(S - wave, Deuteron) = vJT = 0, S = 1) - 8u t(T = 0, S = l) 2 ^ 

Den 

where the energy denominator was adjusted to Den=338 MeV to ensure the correct binding energy of the deuteron. 

The precision of variational estimates for E ( used in the next sections) may be checked against numerical solutions 
of the Schrodinger equation. It is about 0.3 MeV, compared with the overall binding of ~ 50 MeV. The width of the 
state is calculated as 

T/2 =< X nn \ImV K \ Xnn > ■ (26) 



B. Interactions in the decay channels 



The decay channel Y,ir coupled to the basic KN channel is now introduced explicitly. The wave function at each 
scattering center has two components one in the KN the other in the £-7r channel. The scattering amplitudes are two 
dimensional vectors ipi — ► [i/jf^ , ipf] at each nucleon. Multiple scattering equations given in the previous section are 
now changed accordingly. One has 

+ t K ' K G K K + t K " G 7 ^ ij>Z = (27) 



ipl + G 7 ^ V2 + 1*' G V# = (28) 

and an analogous pair with 1 <-> 2. The notation has been changed to describe channel indices and the 2x2 scattering 
matrix /. The determinant related to these equations gives the complex eigenvalue p(xi) in the KN channel. The eigen- 
equation is now more complicated. Introducing a new notation in channel indices U a,b = Q a > a t a ' b the determinant 
becomes 

D =[(1 + U K ^ K )(1 + U 7r ^) -U 7 '' K U K ' 7 <][{l~U K > K ){l-U 7 '' 7 ') -U K ' n U T ' K ]. (29) 
The D = condition is more transparent close to the singularity in the case of a scattering amplitude given by 

; ~ E-E + iT/2- m 

Consistency requires the width to be T/2 = Pui^n) 2 , where is the momentum in the decay channel. The singular 
term p0[) permits one to find a solution of Eq. (|29[) in a fairly simple form. It is presented below in the limit of zero 
range KN ( and E7r ) force. The binding energy 

Re E — E — (7 K ) exp(-p/r) - (7J (31) 

r r 

becomes larger than the binding of the resonance but the collisions in the decay channel indicates oscillations. This 
oscillatory behavior is also seen in the width of the system 

Im E = -(7x) Pit [1 H (Ik) exp(-p/r). (32) 

p w r r 
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The effect of KN scattering represented by the second term enlarges the width as pp is negative. The contribution 
from multiple scattering in the decay channel is sizable in general but it oscillates and may under some conditions 
reduce the total width. That is an effect of interference in the decay channel. Scattering in the decay channel turns 
out to be constructive in the KNN case but it is not necessarily so in some heavier systems. 

C. S and P wave interactions 

With the KN interactions allowed in both S and P waves the scattering equation (fTTj) is a 4 (g) 4 matrix equation 
relating four amplitudes ^/>j. The amplitudes for S waves are now denoted by tpf,^. For P wave interactions the 
corresponding amplitudes are vectors. As there is only one vector in the NN system, the relative separation, the P 
amplitudes are chosen to be Tip p and rip p . The scattering is now described by three types of propagators G a ' related 
to consecutive collisions in the (S, S) (S, P) and (P, P) waves. The scattering equations are 

tpl + f s G ss V2 - f s G sp r 2 $ = (33) 
V>2 + f G ss ipt + f G sp r 2 < = (34) 
ifi + f p G pp tpl + f p G sp V2 = (35) 

1% + f G pp - f p G sp = 0, (36) 

where the propagation in between two P wave interactions is described by G pp = G P q + r 2 G^?. Indices numbering 
the nucleons have been suppressed. The propagator G ss is given in Eq. (fT2|) and explicit formulas for G sp , G P q , G p ^ 
may be found in the appendix. All these functions are regular in the r — > limit. The determinant D of this system 
factorizes into two terms 

D = D S Dp, (37) 

where 

Dg = (1 + G ss f s ){l - G pp f p ) - G sp r 2 f s f p , (38) 

Dp = (1 - G ss / s )(l + G pp f p ) - G sp r 2 f s f p . (39) 
Let us consider the solution of D$ = close to the A(1405) resonance. It is given by an equation analogous to (|T6|) 

r 2 (G sp ) 2 fP 

E = E*-b 2 G ss (r, P (E) [1 - y^^r ]• (40) 

The second term in parentheses describes the effect of P wave interactions. At energies below £(1385) the amplitude 
f p is negative and generates an additional attraction. 

Isospin symmetry simplifies the algebraic structure of the scattering equations which are (see next sections) expressed 
by appropriate isospin combinations of the isospin KN scattering amplitudes /. Equations (|33H36[) allow for a simple 
symmetry of the total KNN wave function. Thus, under the eigenvalue condition (j3"5| the coordinate wave function 
becomes symmetric with respect to the exchange of nucleon coordinates. This condition allows solutions in terms of 
two amplitudes ip s = ipf = -01 an d ip p = ipf = — ^f- Wave functions for the KNN system have the form 

*(r,x) = X NN(r)[G(x - r/2) + G(x + r/2)]^ s + x NN (r) T^ x [G(x - r/2) - G(x + r/2)]^ P) (41) 

where x is the meson coordinate in the NN center of mass system, XAw( r ) is the NN wave function. To make this 
formula more transparent the zero range force limit is taken. The two terms in Eq. (|41[) follow the KN interactions in 
S and P waves. The weight of the P wave contribution is given by 

^ = F (42) 
V> s 1 - f p Gpp ' K ' 
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which becomes dominant close to the zero of the denominator in this equation. At large NN separations it happens 
almost at the singularity in f p . In this region the lowest energy, symmetric, solution of Dg = is given essentially 
by the situation 1 — G pp f p w 0. Such a solution exists for r > 1.6/m in the proper quadrant of the complex 
momentum. This implies that, at large separations, it is energetically profitable for the KNN system to exist in the N 
£(1385) configuration, with the nucleon and £(1385) weakly repelling each other. At shorter distances the condition 
1 + G ss f s w determines the attraction generated by A(1405). Despite repulsive effects of the P wave interaction 
such a solution yields the strongest binding, since a large piece of the binding energy is hidden within the structure 
of £(1385). 

The KNN system is built on short range KN and NN correlations and the ^(r, x) contains a large number of partial 
waves coupled to zero total angular momentum. In Jacobi coordinates (r, s) the L r ® L s decomposition of the first 
term of ^(r, x) is mainly S <8> S. Both terms involve even values of L r and the spin-isospin structure of the NN pair 
is either Inn — 0, Snn = 1 or Inn = 1, Snn = 0. Both types of states may be formed. 

Other solutions are determined by Dp = 0. This condition allows amplitudes of different symmetry ip s = ip\ = — 
and ip p — i/'i = "02 ■ I n comparison with Eq. |4lJ the wave functions for the KNN system now become 

#(r,x) = X N N (r)[G(x - r/2) - G(x + r/2)]0 s + X NN(r)^ x [G(x - r/2) + G(x + r/2)]0 p (43) 

are now antisymmetric in the nucleon coordinates and contain odd angular momenta L r ,L s in the L r <g> L s decom- 
position. The spin-isospin structure of the NN pair is either Inn — 0, Saw — or Inn — I, Snn = 1. The NN 
interactions in the Inn = 0, Snn = states are repulsive and do not support any KNN bound states. On the other 
hand, Inn == 1, Snn = 1 states may be formed. 

The results of the two previous subsections may be unified. The notation used in Eq. (|29jl is now extended to 
include the partial wave index in channel KN : U p ' s = G p ' p t s , U s ' p = G s < s t p , U p ' p = G p > p t p ,U s ' s = U K ' K . The 
last equivalence indicates that the P-wave multiple scattering is included only in the basic KN channel. For the 
determinant of scattering equations one obtains D = D$ Dp where now 

D s = [(1 + U K - K ){1 + U*'*) - U*> K U K '"][1 + U p > p ] - [(1 - u^)U p ' s U a ' p ) (44) 

Dp = [(1 - U K - K ){1 - U 7 *- 7 *) - U 7 '^ K U K - 7! ][l + U p - P ] - [(1 + U^)U p ' s U s ' p ]. (45) 

The solutions of the corresponding eigen- value equations retain the symmetries indicated in the previous section. 

Figure 1 indicates typical contracting potentials obtained with an A.Martin solution (see next section and Ref. [l||). 
Asymptotically these reproduce the separation energies of the K meson hidden either in the 5-wave A(1405) or in the 
P-wave £(1385). The real binding is generated by the difference Vk{t) — Vr-(oo) which jointly with Vnn determines 
the nucleon wave function xnn- 

III. KN INTERACTIONS 

The coupled multichannel KN, £71", Air system is the easiest to describe in terms of the K matrix related to the 
scattering matrix T by the algebraic Heitler equation 

f = K- KiQT, (46) 

where Q is a diagonal matrix of channel momenta in the CM. system. Early parameterizations involved constant 
if -matrix elements chosen to fit the scattering data. Later these were improved by an effective range expansion. As 
the data were (and still are) poor such fits were supplemented by additional consistency conditions formulated in 
terms of dispersion relations [lH, [2(| ■ Such solutions can be tested above the KN threshold and to some extent in the 
£-7r channel. For the dominant isospin interactions there are two types of solutions. These are given in Table [J in 
terms of the inverse M = K^ 1 matrix which in turn determines the scattering matrix 

f- 1 =M + iQ. (47) 

Extrapolations into the complex energy plane display a similar A(1405) pole position. However, the physics in both 
solutions indicates different interplay of the main KN with the hyperon pion channels. The position of the singularity 
is given essentially by the attractive and, in both cases, large Kkn,kn element. This allows one to interpret A(1405) 
as a KN quasi-bound state. In principle there exists an alternative possibility - the A as a quark state. If this is the 
case it may be introduced into the K matrix as an external pole in K ~ 1/(E — E*). However, the scattering data 
exclude such a term or limit it to a very small contribution [19| . Amplitudes below the KN threshold may be tested 
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FIG. 1: Contracting potential Vk(t) in the NN, 7jvjv = 1,Iknn = 1/2 state, for the symmetric meson wave solutions. The Vs 
line shows Re Vk(t) for S wave interactions. The Vsp line shows Re VkM for S + P wave interactions. The upper curve — F 
shows 2 /m Vk(t). These results are based on the A. Martin amplitudes [la ]. 

indirectly, either in the elastic T,tt channel or in the KN — ► £7r transitions on bound nucleons (28|. These reactions 
support the bound state interpretation but are not very restrictive on the position of the singularity. In particular, 
the analysis of Dalitz and Deloff [22j shows that several models offer comparable descriptions of the T,n data in the 
resonance reg ion. The M matrix model given in the DD column of Table U is only slightly favored by the authors of 
ref. [2^| . (3ll |! The KWW column in Table U comes from a quasi-relativistic separable potential model. It belongs to a 
second type of solution and was based on the B. Martin - Sakitt solution. In this work we use off-shell extrapolations 
of both types of solutions. 

TABLE I: The semi-phenomenological 1 = KN scattering parameters. First three lines give M matrices at the KN threshold 
[fm -1 ]. Next two lines give the A(1405) pole position (E* , T/2) in the complex energy plane [MeV]. The KN scattering length 
a + ib and amplitudes at 100 MeV below the threshold Re F_ioo + Urn f-ioo are given in units of fm . The last column 
KWW* corresponds to the KWW model modified to change the A(1405) parameters. The first solution AM, will be referred 
to as type one, and the later solutions BM, DD, KWW and KWW* as type two. 



solution 


AM [19] 


BM[20J DD[22] KWW[11] 


KWW* 


Mkn.kn 


-0.07 


-1.21 


-1.136 


-1.27 


-1.27 


MkN,tvT. 


-1.02 


1.53 


1.254 


1.50 


1.50 




1.94 


-3.05 


-2.205 


-3.05 


-3.05 


E* 


1411 


1415 


1404.9 


1409 


1405 


T/2 


17 


13 


26.6 


22 


24 


a a 


-1.70 


-1.55 


-1.54 


-1.52 


-1.52 


b 


0.68 


0.58 


0.74 


0.60 


0.60 


Re F-ioo 


-2.34 


-1.33 


-8.132 


-2.64 


-5.38 


Ira F_ioo 


0.10 


0.06 


1.11 


0.22 


0.86 



A. The separable off-shell extension 

The three-channel or two-channel separable model is used here to extend the phenomenological S-wave KN interac- 
tions off the energy shell. This method is standard in momentum space but here we have already used the coordinate 
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representation. In a single KN channel case the potential equivalent to those of Eq. ([3]) is described by 

V(k,k') = Xv(k)v(k'). (48) 

The Yamaguchi form-factors v(k) — k 2 /(k 2 + k 2 ) are convenient to perform explicit analytic calculations in both 
representations. The related Fourier transforms are of Yukawa form v(r) — k 2 exp(— «x)/(47rr) and are normalized to 
delta functions in the limit of zero range forces. The off-shell scattering amplitude becomes 

_ v(k)v(k') 

fKN- X - 1+G{E y (49) 
With the non-relativistic form of the kinetic energy E^m = Q 2 /(^^kn) one obtains 

G{E) = f dr J ( 2 t)2 = (50) 
J 27H(t j — q J ) 2(1 — iq/«;)' s 

At the KN threshold energy E t = Mk + Mjv, the standard convention requires one to define the scattering length as 
a + ib = — fKN{Et)- Below the threshold 1/A + G is forced to have a zero corresponding to the A(1405) state. The 
complex momentum at this point is 

Pb = [^fV 2 - «• (51) 

The next step in order to improve the absorptive part is to guarantee that it vanishes below the £7r threshold. That 
is achieved by scaling the absorption strength Im A by a phase space factor f g = q^(E) / q^i^Et) where q^ is 
the momentum in the decay channel. The values As = — 0.602 exp(i 0.12 f ) fm and k — 4.5 fm^ 1 give a good 
reproduction of the PDG recommended E = 1405 and r = 50 MeV values (22l.l25j|. 

This one channel amplitude may serve as a guide, but to describe finer details and for a better comparison with 
the scattering data one needs multi-channel separable models. Below, two types of multi-channel reaction matrices 
are extrapolated off the energy-shell. 

o One solution has been given by Krzyzanowski et al. [ll| in terms of a quasi-relativistic multi-channel separable 
potential. The G{E) used there differs from the solution (|50|) by an invariant momentum phase space and the use 
of quasi- relativistic intermediate meson energies E = q 2 /(2M N ) + ^m 2 +q 2 . This solution was motivated by the 
early B. Martin and Sakitt [2C| K- matrix (BM in Table j]) and as may be seen in this table it offers similar on-shcll 
parameters. 

o Another solution is based on the commonly used A. Martin's K-Matrix [19(. However, in the decay chan- 
nels this matrix is not well reproduced by simple rank one separable potentials. Instead, we use the extrapolation 
K°ff(k,E,k') = v(k) K on v(k'), with the Yamaguchi form-factors and k — 4.5 fm -1 as obtained above in the one 
channel case. 



B. P-wave KN interactions 



To account for the P-wave interactions dominated by £(1385) the scattering amplitude of equation © is generalized 

to 

fKN = fs + fp = fs + 2kk'/j+ (52) 

The last term is a consequence of the j = I + 1/2 total spin of £(1385). It involves an I + 1 factor instead of 21 + 1 
typical for spin zero situations. The omitted piece contains spin flip amplitudes and is expected to be small in the 
few body context. The term is described here by a separable single-channel K - matrix which in the KN channel 
is given by 



where the form-factor is 



« 4 
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and E a is a phenomenological parameter which determines the position of £(1385). The width T n is strongly energy 
dependent 

= 7L<7L + 7L?L- (55) 

In these equations the g's are the channel momenta and 7 are couplings of the resonance to the T,tt and Att channels. 
The latter are determined by the experimental decay width of about 36(±2) MeV and the S7r branching ratio of 
0.13(± 0.01 ) For the coupling to the KN channel the SU(3) value 7^/7^ = 2/3 is taken. This value is 

consistent with the experimental result of Brown, which yields 0.57(±0.18) [23j |. 
The coupling to the KN channel generates the off-shell scattering amplitude 

f P (k,E,k') = 2kkV(fc)Mfc , bLv =kk'v P (k)f P (E)v P (k'). (56) 

Let us notice that below the if AT threshold the integral in Eq. (p)6")) deforms significantly the shape of the resonance 
profile. The range parameter k = 3.8 fm _1 and E a = 1505.2 MeV are used to reproduces the profile tested experi- 
mentally by Cameron et al. 24] in the An channel. For further applications the coordinate representation is needed, 
which is given by equation 

/(x, E, x') = A P *d v P (x)fp(E)v P (x')W. (57) 
in terms of the Fourier transforms of form- factors (|54[) . 



IV. FEW NUCLEON SYSTEMS 



The procedure presented in the KNN section is now extended to systems consisting of several nucleons. Practical 
calculations are done for three and four nucleons. At first the multiple scattering equations similar to Eqs. ([3] -[5]) in 
the previous section are solved in fixed nucleon systems. The bound K-meson wave function \k is a solution of 

Xif(x,xi....x„) = -SiS/3 / dy — - — : : — v(x-x i ) Q A Q/3 v(y - Xj)^ Xk(y, Xi....x n ), (58) 

J 47r | x - y | 

where indices a, (3 denote channels and partial waves of the meson-baryon pair. An index related to the symmetry of 
X is suppressed. By analogy with Eqs. (6) and (7), equation (|58l) may be reduced to a matrix equation for the wave 
amplitudes defined at each scatterer as 

i>f = E/3 J dy \ afi v(y - xO/3 x K (y, x ; ). (59) 

The kernel of the scattering equation can be now expressed in terms of scattering amplitudes /" ,/3 at each nucleon i 
and propagators describing the passage from nucleon i to another nucleon j. The latter are given by 

n«>P [a a l \ ex P( ik I x -y I) / \ 

G it 'f = J dy dx v(x-Xi) a ^ ^ — v{y-x. j )p. (60) 

The procedure explained in Sect. II leads to a set of linear equations 

C + S^^/^G^VJ^O, (61) 

which are solved numerically. This matrix equation is simplified as the G's are diagonal in channel indices and the 
/ are diagonal in the partial wave index. Still, the corresponding determinants are complicated algebraic expressions 
involving functions G and /. Numerical solutions become a difficult problem. It has been solved in the following, 
approximate way. The determinant consists of many terms which are arranged according to the number of collisions. 
With up to four collisions in the main channel we retain the structure found in the KNN situation and the determinant 
D of this system is presented as 

D = 1 + Y. pairs (D S Dp — 1) + Ohigher orders, (62) 

where D$ and Dp are defined in Eqs. (|4"4"]l and (14"5"1) . The main term is composed of collisions in the KNN subsystems 
which allows one to keep track of the wave function symmetry. The terms of higher order in / are dropped. 
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The solution of the full K-few-N bound state problem is given by the equation 

i i i,j 

Again we assume the wave function to be given by eqs.([2]) and ([55]) i.e. Xk(%i x\..x n )xN{xi, .., x n )). Projecting 
Eq. (|63[) on tpK one obtains the Schrodinger equation for the few nucleon wave function 



XN(xi) = j dx XK(x,Xi)^(x,Xi) (64) 

in the form 

[E - E c ( Xi ) + ^ ~ v nn]xn + AE kmXN = 0, (65) 

i 

where AEkin is a correction to the nucleon kinetic energies. As in Eq. |22|) it is given by 

AE kinXN = J dxd^Kdi* (66) 

and as before turns out to be very small due to angular averaging and sign changes in both the derivatives. As 
discussed in the KNN situation this correction has been dropped. In deriving equation (f65|) a special form (|58[) of the 
meson wave function is used. As in Eq. (|21|) it satisfies the relation 



A xXK = Y. A ** K - ( 67 ) 



In the next step, equation (|65p is solved by a variational method with the NN potential of Argonne [l5| . The trial 
wave function is of the form 

XN= J] t 1 cx P(- A c(*i - xj) 2 )] exp(-A, | x; - xj |)/ | x, - xj |, (68) 

where A c , A; are variational parameters. This form is chosen to give the correct asymptotic limit for large | Xj — xj | 
and also gives a vanishing wave function for small | x ; — xj | as expected for a strong repulsion in the NN potential. 

V. RESULTS 
A. KNN 

In this section the calculations of KNN levels are presented. The sensitivity to KN input parameters is studied and 
the states of different symmetry are discussed. 

Contracting potentials Vk (r) were calculated with several solutions for the phenomenological S-wave KN reaction 
matrices presented in Table |TJ The solutions of second type may be well fitted with a rank one separable potential. 
Here, the calculations are done with the quasi-relativistic model of This model is based on the if-matrix of B. 
Martin [2(| . Numerically it is fairly close to the separable potential of Ref . [2l| . For the the first type of solutions in 
Table U- due to A. Martin [l9| no satisfactory rank one separable approximation is found. This difficulty is related 
to the large effective range parameters involved in this K matrix. In order to retain the physics involved, a simple 
off-shell extension is adopted Ki j — > v(ki)Ki t jV(kj). The Yamaguchi form-factors have been used and the inverse 
range parameter k was varied over the range of 3 — 6 fm . The actual value of k affects the multiple scattering 
via propagator G ss (k,r) of Eq. (fT2"]) . Larger values reduce the form-factor v(k) but enhance the significance of the 
small r region in G ss (k,r). On average these two effects balance very well and one finds a very weak (~ 1 MeV) 
dependence of the total binding energy on the actual value of k. The results given in Fig.l and in the tables which 
follow are obtained with k = 4.5 fm -1 . 

The energies of the most strongly bound KNN, I tot = 1/2, Inn = 1, quasi-bound states are given in Tables HT1 and 
IIIII The first table describes several steps of the approximation while the second table indicates the dependence of 
binding on the KN input parameters. 
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The first line in Table [IT] is determined essentially by the effects of A(1405) excitations described in the elastic 
channel only. The second line describes additional effects due to multiple scattering in the £7r channel. The other 
two lines include the P wave interactions. The energies of I to t — 1/2, Inn — 1 states given by the S wave interactions 
and described by multiple scattering in the single, KN, channel span the region of 30-50 MeV. These results are 
consistent with the findings of Akaishi and Yamazaki [9( and Dote and Weise [12| obtained with different methods. 
The differences within this range are due to a different KN and/or NN input. As seen in the second and fourth rows, 
significant changes arise with the explicit inclusion of the multiple scattering in the £, n channels. The binding rises 
by 10 to 20 MeV and the effect of collision broadening is large. 

TABLE II: Binding energies and widths [in MeV] of the KNN, I tot — 1/2, Jjvjv — 1 space-symmetric states. The results on the 
left are based on [ly] parameters, the results on the right follow KWW [llj parameters discussed in Table [T] The first column 
specifies the channels explicitly involved in the multiple scattering and meson-nucleon partial waves. Rrms is the radius mean 
squared of the N-N separation [in fm ] . The last line is obtained with the simplest separable potential discussed in the text and 
the Ikn = 1 amplitudes from AM. 



solution 


AM 


19J 


KWW [llj 




E B 


r 


Rrms 


E B 


r 


Rrms 


KN; S 


27 


36 


3.1 


35.5 


37 


2.4 


KN,Ztv; S 


37 


42 


2.5 


43.1 


47 


2.1 


KN; S,P 


49 


36 


3.7 


49.7 


36 


3.3 


KN, £-7r; S,P 


52 


37 


2.9 


56.5 


39 


2.3 


KN; S 


47 


47 


2.3 









Table Hill shows binding energies obtained with the "canonical" A pole position E = 1405 MeV [13, HH which 
is lower than the position obtained in other parameterizations. The result given in the second line of this table 
is comparable to the results obtained, with a similar input, by Schevchenko et al.fljjj]. The latter work employs a 
superior Faddeev technique, but a more detailed comparison of results is not easy since that calculation uses a rank 
one separable potential to describe the NN interactions. 

TABLE III: Bindin g en ergies and widths [in MeV] of the KNN, I to t = 1/2, Jjvjv = 1 space-symmetric states. These results 
are based on KWW 41] parameters modified to set the pole of A(1405) at 1405 MeV and the width at 48 MeV and given in 
Table |T] (KWW*). The first column specifies the channels explicitly involved in the multiple scattering and the meson-nucleon 
partial waves. Rrms is the radius mean squared of the N-N separation [in fm]. 



solution KWW* 






E b 


r 


Rrms 


KN; S 


50 


51 


2.05 


K7V,Ett; 5" 


71 


85 


1.81 


KN; 8,P 


65 


43 


2.09 


KN,T,it; S,P 


78 


60 


1.88 



The inclusion of resonant P wave interactions increases the binding by some 10 MeV. There is some room for 
different values as the experimental K7V£(1385) coupling is not certain. However, the main effect of £(1385) is a 
change of structure in the KNN systems. A sizable portion of the binding energy is contained in the structure of this 
resonance. On the other hand the system is dissolved as the inclusion of P waves enlarges the NN separation and the 
formation of £(1385) is essentially a peripheral effect. The KN correlations for r > 1.6 fm are mostly of the £(1385) 
type. The other effect of P wave interactions is a formation of additional KNN states. These are given in Tables HVl VI 
and discussed below. 

The energies of KNN quasi-bound states with Inn = given in Table HVl are determined essentially by the £(1385) 
excitations. Let us notice that the result is unstable against the KN input. The state is still more likely to exist with 
a lower value of the A(1405) energy. In any case it is a very loose structure that might be a quasi-bound or a virtual 
state. 

The energy of an asymmetric quasi-bound state is given in Table fVl It is determined essentially by the £(1385) 
excitations. The table in Appendix B indicates that the K~nn state has the largest possible £ component which 
offers the strongest £(1385) — N attractive potential. It reaches a maximum depth of about —10 MeV at 1 fm 
distance, but it is not strong enough to overcome the NN P-wave barrier and generate a quasi-bound state. To obtain 
real binding, assistance from the Ikn — 1 S- wave state and the NN P-wave attraction is necessary. Thus, the NN 
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TABLE IV: Binding energy and width [in MeV] of the KNN, Itot = 1/2, Inn = space-symmetric states. Results on the 
left are based on A.Martin parameters and do not support a bound state (no b.s.). The results on the right follow the KWW 
parameters. This result involves S + P wave interactions and external interactions in two meson-nucleon channels. R r ma is the 
radius mean squared of the N-N separation [in fm]. 



E B 


r 


Rrms 


E B 


r 


Rrms 


no b.s. 






47.1 


36 


~7 



TABLE V: The binding energy and width [in MeV] of the KNN I tot = 3/2, Inn = 1 space-asymmetric states. In the NN 
subsystem 2S+1 Lj — :i Pi. The last column gives the radius mean squared of the N-N separation, [in fm]. 





Innk 


Inn 


E B [MeV] 


T[MeV] 


-Rrms [fm] 


K nn 


3/2 


1 


48.5 


36 


4.9 



interactions repulsive at large distances in the Inn = 0, Snn — 0,Lnn = 1 waves do not support bound states. 
However, such states may be generated by Inn — 1 , Snn = 1,Lnn = 1 interactions. Here, the analysis becomes 
more subtle as the NN interaction is strongly spin dependent. The energy given in Table fVl corresponds to J = 2 
( 3 P2) wave in the NN subsystem where the interaction is the most attractive. This KNN system is large and loosely 
bound, by about 1.5 MeV in the £(1385) — N configuration. In this calculation the S wave parameters come from 
AM [191 ] and the calculated energy is uncertain, as the involved K~n parameters are poorly known. The experimental 
detection would be very difficult, nevertheless, a more precise analysis of the spin and isospin structure of such states 
is of interest in the context of K~ D atoms and will be performed elsewhere. 

B. KNNN and KNNNN systems 

The discussion of these systems is limited to the states of the simplest symmetry. The fixed nucleon model generates 
a contracting potential which in KNNN systems may be, to a good approximation, presented in the form 

E c {R x ,R y ,R z ) = -V N NN[l-Ccxp(-\ s {R x + R y + R z ))}[exp(-\ l R x ) + exp{-\ l + R y ) + cxp(-\ l R z )}-V{™), (69) 

where R x , R y , R z are the inter-nucleon distances. The short range behavior at the triple coincidence may be obtained 
analytically and C = 0.42, other parameters being numerical. 

With the KN parameters of Refs.[llj],[lll an d Itot = the parameters are obtained in the range Vnnn ~ 150 — 200 
MeV, A s ~ 4.5 fm , A/ - 1.8 - 1.9 fm. For I tot = 1 one has V NNN - 50 - 60. V(oo) is the binding of KN into A(1405) 
in the S wave case or £(1385) in the S+P wave case. For I to t = the corresponding binding energies are given on the 
left side of Table fVTl These numbers may be compared to the simplest version of this model - the S wave interactions 
described by the single KN channel - which produce 91 MeV binding. 

The modified version of the KWW model with parameters from Table Q] fixed to set the A(1405) energy to 1405 
MeV yields much stronger contracting forces, Vnnn ~ 250 — 350 MeV and A; ~ 2.1 — 2.3 fm. The states indicated on 
the right side in Table IVT1 are bound very deeply. The basic NNN systems obtained with our variational wave function 
are over-bound by about 2 MeV and this value has already been subtracted from the numerical KNNN energies. 

TABLE VI: Binding energies and widths [in MeV] of the Itot = 0, KNNN, space-symmetric states obtained with the two channel 
KN , E-7T channel multiple scattering formulation. The results on the left are based on the modified KWW [ll[ parameters with 
A(1405) energy set to 1405 MeV. The first column specifies meson-nucleon partial waves involved. The widths do not describe 
non-mesic capture modes. 





E B 


r 


E B 


L 


s 


103 


29 


142 


25 


S + P 


119 


23 


153 


21 



There may exist a number of states in the KNNNN systems. In Table IVHll one finds only the states with the simplest 
symmetry, which involve wave functions symmetric under exchange of the nucleon coordinates. In the absence of the 
K meson the basic a particle structure is used and only the S wave NN interactions are included. With the tensor 
interactions described by Eq. (j2"5"|) this a system is over-bound by about 10 MeV, and this value has been subtracted 
from the calculated KNNNN levels. 
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TABLE VII: Binding energies and widths [in MeV] of the hot = 1, KNNN, space-symmetric states obtained with the two 
channel KN , E-7T channel multiple scattering formulation. These states are formed as a result of the P wave interactions with 
some assistance of the S wave attraction. 





E B 


r 


S + P 


63 


38 



TABLE VIII: Binding energies and widths [in MeV] of the KNNNN , space-symmetric, Stot = 0, hot = 1/2 states. See captions 
to Table IVll 
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E B 
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121 


25 


170 


10 


S + P 


136 


20 


172 


10 



C. Level widths 

Level widths are calculated as twice the expectation value of Im Vk ■ The KN resonant lifetimes are strongly energy 
dependent, being very short at the resonance, becoming longer below the resonant energies and staying infinite below 
the thresholds of the decay channels. This trend is reflected by Im Vk in Fig. [T] The energy dependence contained 
in the amplitude Jkn{— Eb — E rcco u) of Eq. (P) is traded into the space dependence /kn^r)- These two types of 
averaging give fairly close results provided the final binding energy is located well above the threshold of the decay 
channels. Let us indicate some consequences of this relation. 

The states generated by the P wave interactions given in Tables ITV I V[ and I VIII correspond to a fairly loosely bound 
£(1385) and the widths of quasi-bound states are essentially equal to the width of the £(1385). This comes as a result 
of the peripheral binding and weak effects of the collision broadening in the P wave resonances. In these states, the 
Vk underestimates slightly the average — Eb — E rcco i\ and the real widths might be smaller. For the binding energies 
in the range 60-90 MeV the Vk is too small at large distances and too large at small distances with a reasonably good 
average. Let us notice that the level widths generated by the S + P interactions are smaller than those generated by 
the S waves alone. This is due to three factors: the width of A(1405) is larger than the width of £(1385), the collision 
broadening in P waves is small and the systems due to the S + P interactions are less compressed. Let us also notice 
very strong sensitivity to the input KN amplitudes. The few examples of Im -Fioo given in Table U and the differences 
of the widths in the I to t — 1/2, Inn — 0, KNN state visualize this point. 

In cases of very large binding, in the range of 120-200 MeV, one (£7r) or two (£7T, A7r) decay channels are blocked 
and the widths calculated here are over estimated. Such a situation is likely to happen in the K-a state. To account 
for that effect, the calculation of the contracting potential was repeated in an optical potential manner. So, the 
momenta in the decay channels £7T, A7r were related to the binding energy Eb and allowed no outgoing waves. Such 
a calculation results in a stronger binding. In the KWW* model one obtains binding of 220 instead of the 170 given 
in Table IVIIH The real decay width is now given by the multi- nucleon capture mode. 

The multi-nucleon captures are initiated by the non-mesic KNN — > Y N mode and the branching ratio for this 
process is known from emulsion studies to be about 20% in light nuclei [30| . The emulsion data are obtained with 
stopped K mesons and pertain to the nuclear surfaces. An extrapolation in terms of a characteristic nuclear densities 
p and two body phase space L 

Tmulti — L g 2 7 (70) 

for this decay was attempted in Ref. A constant 7 may be fixed to the emulsion branching ratio and a 20 MeV 
level width in the nuclear matter at 90 MeV binding was obtained. In the strongly bound, few-body systems the 
kinematics of the decay is different since the residual nucleons also take sizable recoil energies. Roughly, for a three 
body decay L ~ Q 2 where Q is the decay energy. Again, an extrapolation from the emulsion data in terms of the 
available phase space and the involved nuclear density yields non-mesic capture widths in K-a in the 10-30 MeV range. 
These estimates are somewhat larger than the 12 MeV obtained for the KNNN system in Ref.Q. Unfortunately such 
extrapolations are uncertain as the energy dependence in 7 might be large and the Q value is not known. Help from 
new experiments is necessary to settle these questions. 
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VI. CONCLUSIONS 



In this paper a new method to calculate the deeply bound KNN, KNNN and KNNNN states has been presented. 
The calculation consists of two steps. First a wave function involving strongly correlated K-N subsystems is found in 
a fixed nucleon approximation. This step also allows one to find potentials due to the K meson which tend to contract 
the inter-nucleon distances. Next, these correlated wave functions and contracting potentials are used as input in the 
Schrodinger or variational calculations for the K-few-nucleon binding. 

The lowest binding energies based on a phenomenological KN input fall into the 40-80 MeV range for KNN, 90-150 
MeV for KNNN and 120-220 MeV for Ka systems. The uncertainties are due to unknown KN interactions in the 
distant subthreshold energy region. 

We obtain at least partial answers to the basic questions presented in the introduction. 

(a) The binding mechanism: the dominant mechanism of the attraction is related to the A(1405) state. This fact has 
been known for a long time. In addition, it is found here, that the £(1385) contributes significantly to the structure 
of the K-few-N bound states but much less to the actual binding energies. The bound states are built from the 
strongly correlated KN subsystems. At central densities these correlations resemble the A(1405) and at peripheries 
the correlations are made by the quasi-free £(1385). Sizable fractions of the binding energies are contained in the 
KN correlations. One consequence is that even with the strong bindings the nucleon densities are not dramatically 
enhanced as in Ref. 0] but can become a factor 2-4 larger than the standard nuclear matter density p . 

The presence of £(1385) resonances in the few nucleon systems generates new states. These are predominantly 
P-wave states or states built upon the P-wave N-N interactions and are usually broad and loosely bound developing 
long tails built from the £(1385) correlation. 

(6) The control of technical questions: the choice of correlated wave functions removes the difficulties related to 
the uncertain K-N interaction range and allows one to use realistic N-N interactions. The recoil energy of the KN 
subsystems with respect to the residual nucleons is described only in an average sense. This seems to be the weakest 
part of this method. 

(c) The widths are related to the lifetimes of the A(1405) and £(1385) enhanced by the collision broadenings. Under 
the phenomenological KN interactions ( the A pole located at ~ 1412 — 17 i MeV ) the K-few-N states are ~ 40 MeV 
wide. 

On the other hand, the models with the A pole located at ~ 1405 — 25i MeV generate K-few-N states which are 
more deeply bound. These may be either very broad or quite narrow. With the binding energies in the 60-80 MeV 
range ( the KNN case) very broad - up to 90 MeV - states are obtained. However, with the bindings of 120 MeV 
(KNNN) or 170 MeV (KNNNN) the single nucleon decay modes are effectively blocked. The widths are strongly 
reduced and the main decay modes are due to multi- nucleon K captures. These widths are hard to predict, a simple 
model suggested here generate widths of about 20 MeV. 

A simple physical picture emerges from this approach. The mesons are strongly correlated to slowly moving 
nucleons. The correlations are of the A(1405) type at large densities, and of the £(1385) type in the peripheries. Each 
K,N pair has a good chance to stay also in the £,7r form. The structure is rather loose as sizable fractions of the 
binding energies are hidden in the short ranged correlations. 
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APPENDIX A: PROPAGATORS 

Several formulas for kernels of multiple scattering equations are collected in this appendix. 

For Yamaguchi form-factors, the propagators Gij yield analytic expressions. Thus for two consecutive S-wave 
interactions one has 

i r k 2 + k 2 i 

G(r,k) ss = v(k) 2 expUkr) — exp(«r) — r— exp(— nr) , (Al) 

47rr L 2k J 

wherer = Xi — Xj, and the indices i, j referring to the nucleon sites are suppressed. For an initial S wave scattering 
followed by a P wave scattering G becomes a vector 

G(r,k) sp = rG sp (r,k), (A2) 



1 



r 2 {K 2 + k 2 ) r 3 (K 2 + k 2 ^ 2 ^ 



G sp = - — v(kY exp(ikr)(ikr - 1) - cxp(«x)[l + nr H — H ^— — ] . (A3) 

For two consecutive P wave interactions the propagator is a tensor of the form 

G(r, k) pp \ nm = v{k) 2 [5 nm G P J + r n r m G P T P ]. (A4) 
These functions may be expressed in terms of basic integrals 

47T f K 2 1 



which give by recurrence 



exp(ifcr) - exp(«:r) k 2 + k 2 
3 i (r) = , g 2 (r) = gi (r) ^— exp(-Kr) (A6) 



gz{r)=gi{r) ^ (1 + Kr)exp(-Kr), g 4 (r) = g 3 (r) — (1 + nr H — )exp(-Kr) (A7) 



and finally 



rPP _ 9i{r)' pp _ gj{r)" _ g A {r)' 



(A8) 
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APPENDIX B: ISOSPIN SYMMETRY 



It is assumed here that the isospin is conserved in the quasi-bound states of K mesons. In the lowest S-wave states 
of the KNN systems the isospin wave functions may be built upon iso-singlet or iso-triplet NN states. From the 
experimental point of view, the most interesting one seems to be 

^{ /2 = {{NNyK} 1 ! 2 = V3/2{{NK}°N} 1 / 2 + 1/2{{NK} 1 N} 1 / 2 , (Bl) 

1/2 

where in 1 S r 1 / the upper index denotes isospin I nuc i and the lower index denotes the spin of the two nucleons. On 
the right side the isospin content in the KN subsystem is given. This state is a mixture of K~pp and K°np and is 
frequently named K~pp since it can be experimentally accessed via this entrance channel. The NN spin in this state 
is S — and the effective KN interaction amplitude obtained from Eq. (|Bl[) becomes 

fKN = 3/4 f KN + 1/4/V (B2) 

Another KNN state of interest is built upon the NN iso-singlet 

q>l /2 = {{NN}°K}^ 2 = -1/2{{NK}°N} 1 / 2 + VS/2{{NK} 1 N} 1 ^ 2 . 

This state is a mixture of K~np and K°nn which might be reached by the K~np entrance channel 
spin is S = 1 and the effective KN interaction amplitude obtained from Eq. IB3I becomes 

fKN = l/4f°KN + 3/4fk N . (B4) 

The S - wave KN interaction in the \&j|/ 2 state is much less attractive than in the ^>\^ 2 state since the A(1405) 
contribution is reduced. However, this is offset by the strong short range attraction in the NN system due to the 
tensor force. An additional attractive force is due to a larger contribution from the £(1385) resonance. 
Finally one may have total isospin 3/2 states of the type K~nn or K°pp 

m\ /2 = {{JVJV} 1 ^}^ 2 = {{NR^N} 3 / 2 . (B5) 

Those states, involve weakly attractive and uncertain, S - wave KN 1=1 amplitudes. A deeper state can in principle 
be built upon the stronger P - wave interactions. Its existence and the chances for detection present a situation that 
is more difficult than the other cases. 

For the three nucleoli problem we retain the dominant structure of the triton and helium isospin symmetry. The 
KNNN wave function is assumed to be of the form 

* T = -^{{{{AiV} ' 1 ^} 1 / 2 + {{NNy^N} 1 / 2 }^, (B6) 
v2 

where the pair of indices denote spin and isospin of the NN pair. 

Re-coupling to the KN system leads in the total T=0 state to the relation 



(B3) 

Now the NN 



*° = ^lj2{{NNY{NKY + {NN}°{NK}°}° (B7) 
and in this case the KN interaction amplitude is 

= l/2f° KN + l/2f KN . (B8) 

Likewise for the total isospin 1 system 

/ s = l/6/kv + 5/6/^ (B9) 



These amplitudes are collected into the table. 
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TABLE IX: Isospin composition of Kaon nucleon scattering amplitudes. I tot = total isospin, I nuc i = isospin of nucleons, /, = 
KN amplitudes of isospin, i. 



System 


hot 


Inuci 


/if N 


KNN 


3 
2 


1 


h 


KNN 


1 
2 


1 


|/o + \ fl 


KNN 


1 
2 





4/0 + 


KNNN 





1 

2 


5/0 + 


KNNN 


1 


1 

2 


^/o + I/i 


KNNNN 


1 
2 





3/0 + 3/1 



APPENDIX C: THREE NUCLEONS, S-WAVE INTERACTIONS 

The energy eigenvalue is obtained by the simultaneous solution of three equations 

4>i + G{%f'ih + Gi%f'M = (CI) 

r 2 + G' 2 %f'P 3 + G'frf'M = (C2) 

^ + G&rri + G s 3 %rr 2 - 0, (cs) 

which require the eigenvalue condition 

lh. I (/*) 2 [G$G;? 2 + G{%Gi% + G S 3 %G S 3 %] + 2{f') 3 G' 1 %G% 3 G% 1 = 0. (C4) 

This equation is to be solved numerically. A helpful guide to find the symmetry of two physically meaningful solutions 
is the situation of two equal NN separations m = r\ 3 . Dropping the upper indices one obtains the factorized form 

D 3s = (1 - / Gi, 2 )(l + / Gi, 2 - 2 f 2 Gl 3 ). (C5) 

The first factor corresponds to an antisymmetric solution with the meson sticking to two nucleons only. The sec- 
ond factor generates a solution symmetric with the interchange of nucleons 1 and 2. These solutions are a direct 
continuation of the two solutions obtained in the KNN systems. 



